Single atom-scale diamond defect allows large Aharonov-Casher phase 
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We propose an experiment that would produce and measure a large Aharonov-Casher (A-C) phase 
in a solid-state system under macroscopic motion. A diamond crystal is mounted on a spinning disk 
in the presence of a uniform electric field. Internal magnetic states of a single NV defect, replacing 
interferometer trajectories, are coherently controlled by microwave pulses. The A-C phase shift 
is manifested as a relative phase, of up to 17 radians, between components of a superposition of 
magnetic substates, which is two orders of magnitude larger than that measured in any other atom- 
scale quantum system. 
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One of the fundamental predictions of quantum me- 
chanics is the existence of topological phases. A famous 
case still relatively unexplored is the Aharonov-Casher 
(A-C) effect: if a particle with a magnetic moment is 
taken about a line of charge, it acquires a phase which 
is independent of its velocity and only topologically de- 
pendent on its path [lj. The effect is important for spin 
currents in mesoscopic rings 0, Q , vortices in supercon- 
ductors [1, [B| , and it has potential applications in topo- 
logical quantum computing @, Q. The A-C phase, a 
relativistic effect, is much harder to produce than the 
Aharonov-Bohm effect [8j and other better-studied topo- 



logical phases. Experiments using atom beams [3, LKJ, 111 | 
and neutron interferometers [l2| have successfully mea- 
sured the A-C effect but only small phases have been 
produced, of order 0.15 rad or less. In contrast to these 
atom-scale systems, large A-C phases can arise through 
collective effects in mesoscopic rings Here, we pro- 
pose a method to generate and measure a large A-C phase 
through quantum control of a single atomic-scale system: 
the negatively charged diamond nitrogen-vacancy (NV) 
center. 

The NV center is a naturally occurring defect in di- 
amond which has a spin triplet electronic ground state 
with excellent coherence properties, making it a widely 
considered candidate for a quantum bit [13l [3]. In 
our proposed experiment, an NV center, coherently con- 
trolled by microwave pulses, acquires a relative A-C 
phase between its magnetic sublevels as it moves be- 
tween two charged plates (Fig. [Ha)). The experiment 
is made possible by recent developments in coherent con- 
trol of NV centers, which have shown that an NV center 
can act as a highly sensitive magnetometer, detecting 
fields down to nTHz -1 / 2 sensitivity, at room tempera- 



ture [15J, US E3, EI 0, |2fJ. The experiment would be a 
powerful test of the A-C effect, generating an accurately 
measurable A-C phase two orders of magnitude larger 





Figure 1: Proposed experimental setup, (a) Geometry of ex- 
periment; a diamond crystal on a spinning disk between two 
charged plates with a uniform, static magnetic field in the 
z-direction. (b) Energy level diagram of the NV center. 



than those of previous experiments, and could open up a 
new field of applications for coherently controlled solid- 
state quantum systems. 

A-C phase. — Aharonov and Casher considered a neu- 
tral, non-relativistic, spin-^ particle, with magnetic 
dipole moment /i<r, moving in a chargeless region of non- 
vanishing electric field £(x). The Hamiltonian for such 
a particle can be derived from the Dirac Lagrangian in 
the low-energy limit: 
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The A-C effect is usually treated in 2+1 dimensions. 
Some authors use an explicitly two-dimensional wave 
equation 21[ . Others effectively reduce the Hamiltonian 



to two dimensions by requiring the wavefunction to be 
independent of the z-coordinate, and setting £ z — and 
d z £ = [22j]. In either approach, if \tpo) is some solution 
to the Schrodinger equation in the absence of an electric 
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solves the Schrodinger equation in the presence of the 
electric field. The path of integration 7 can be any path 
in the x-y plane that shares the same topology as the 
path of the particle. \ip) differs from |^>o) only by an 
SU(2) phase $(x, y) which is the A-C phase. 

We briefly outline an alternative derivation of the 
A-C effect, which holds for a localized particle of arbi- 
trary spin in three dimensions. For particles of arbitrary 
spin S, the Dirac equation generalizes to the Bargmann- 
Wigner (B-W) equation [21]. Taking the B-W equation 
in the low-energy limit gives the Hamiltonian 
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Where S is the spin operator of the particle, V is a poten- 
tial which may depend on spin, and fj,S is the magnetic 
moment of the particle. We can ignore the last two terms 
of Eq. (j4|) since, due to the pulse sequence discussed be- 
low, they will not contribute to the measured signal. 

Let \ipo;t) be a solution to the Schrodinger equation 
with £ = 0, that describes a highly localized particle 
following a trajectory r(t). A solution to the Schrodinger 
equation with £ =/= can then be written as 



= A(t)U(r(t))\^ ;t), 



(5) 



where A(t) = 1 + X £(r(t)) ■ [x - r(t)],(6) 

he 1 

And U must satisfy an integral equation, 

U(r(t)) = l + ^j\sx £(r(t')) ■ r]U(r(t'))dt' (7) 

whose solution can be expressed as a Dyson series. We 
have assumed that [U(r(t)), V] = 0, and that the particle 
is sufficiently localized that x\ipo;t) w r(t)\ipo;t) and, 
hence, A\ip ) ~ iVjo)- 

The operator U{x) is an SU(2) phase factor, similar 
to the A-C phase factor exp[— i*f>(x, y)] of Eq. (|3|). How- 
ever, whereas only depends on the path's topol- 
ogy, U (x) depends, in general, on the details of the path 
r(t) of the particle. Small movements of the particle 
cause precession of its spin about an axis perpendicular 
to both the electric field and the direction of motion. The 
path-dependence of U arises from the non-commutativity 
of rotations about different axes: as Anandan notes [231 ] . 
we can view S x £ as a non-Abelian gauge field. 

The full SU(2) phase can be explored by tilting the 
plane of the spinning disk about the y-axis. However, 



the non-Abelian effects will be small, as they are of sec- 
ond order in the Dyson series expansion of ([7]). If we 
restrict the particle's motion to the x-y plane, and set 
£ z = and d z £ = as before, we recover the usual, 
path-independent, A-C effect: U(x) = e~ !# . 

A-C phase in diamond. — The NV electronic energy 
level structure is shown in Fig. [Tfb) . We are interested 
in the 3 A ground state, which is connected to the excited 
states by optical excitations. Initialization of the NV 
center is performed by optical pumping into |0), followed 
by Rabi pulses to achieve some coherent superposition 
of the spin sublevels |0) and |1) [26]. The relative pop- 
ulations of the spin sublevels can be observed by optical 
excitation and a measurement of fluorescence. 

Ignoring crystal asymmetries and interactions with nu- 
clear spin, the spin Hamiltonian H s for the 3 A ground 
state can be written as 



H s = D{S Z 



g^sB ■ S, 



(8) 



where D = 2.88 GHz is the zero-field splitting parameter, 
B is the external magnetic field, and g ss 2 is the gyro- 
magnetic ratio of the NV center. Note the NV center's 
negative charge does not affect the spin Hamiltonian and 
hence is decoupled from the A-C dynamics. 

To treat the A-C effect for the NV center, we introduce 
a Hilbert space that is spanned by the position eigen- 
states of the NV center. A general state of the NV center 
can then be written as \ip) = \ips) ® I^r) where \ips) is 
a superposition of |0), |1) and |-1), and \ipR) gives the 
position of the NV center. The total Hamiltonian of 3 A 
state is then 
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where p is the NV center-of-mass momentum, m is the 
NV center's effective mass, and V is a potential which 
acts only on the position space. 

Now consider an NV center moving in the x-y plane, 
in the presence of an electric field £(x), with £ z = and 
d z £ = 0, and a magnetic field, of magnitude B, in the 
^-direction. Let \r(t)) be a solution to the Schrodinger 
equation for the position Hilbert space, such that \r(t)) 
describes the motion of the NV center along some trajec- 
tory, r(t). If the NV center is initialized to 

1-0;* = 0) = (co|0}+ci|l»<8|r(t = 0)) (10) 

where the a are complex coefficients, then after time t 
the state of the NV center will be 

|V>;t) = [c |0) +e -^ Ac{ - r{t)) e- lt ^ D+9B ^/ h c 1 \l)] 
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It is this &ac( x ), the A-C phase in diamond, which our 
proposed experiment would measure. 



We briefly pause to consider the topological nature of 
the AC effect. Boyer 24[ and Casella 25j, among others, 
have contested the appropriateness of the term 'topolog- 
ical'. In particular, Casella noted that the A-C effect can 
be detected without a topological defect and proposed 
the geometry that we use here: a particle in a superpo- 
sition of states passes between two charged plates. (The 
A-C phase in this geometry is often termed the 'Casella 
phase'.) Gorlitz et al. [10| make a strong case for ob- 
serving the A-C phase using this geometry, in which the 
superposed substates are treated as arms of an interfer- 
ometer. The characteristic nondispersiveness and path- 
independence of topological phases remain in any case. 

Experimental proposal. — We envisage an experimental 
setup as illustrated in Fig. [Ha). A diamond crystal is 
embedded at the edge of a disk of radius r that lies in 
the x-y plane and spins, with frequency /, about the 
z-axis. The crystal is oriented so that the NV of inter- 
est is aligned parallel to the z-axis. Two charged plates, 
connected to a constant voltage source, provide a uni- 
form static electric field £ in the x-direction. A pair of 
Helmholtz coils provides a uniform static magnetic field 
B, on the order of Gauss, along the z-direction, which 
lifts the m = ±1 degeneracy jl5| . The experiment con- 
sists of initializing a single NV center into a coherent 
superposition of spin sublevels and observing the rela- 
tive phase that accumulates between the sublevels due 
the A-C effect. As the disk rotates, the NV spin states 
evolve according to . Note that the rate of A-C phase 
accumulation varies sinusoidally in time. Positive (neg- 
ative) A-C phase accrues as the diamond moves in the 
positive (negative) y-direction. 

The NV center is initialized, controlled and measured 
using a 532nm laser, a microwave generator, and a pho- 
todetector, which each focus on particular points of the 
diamond's path as shown in Fig. [2]^a). The laser is used 
for optical pumping and readout and the photodetector 
measures fluorescence. The microwave generator drives 
two coils, located close to the spinning disk, at points A 
and B, which perform controlled Rabi rotations of the 
NV center's spin sublevels. Alternatively, Rabi oscilla- 
tions could be induced optically [27j |. 

The pulse sequence and the corresponding quantum 
state evolution are shown schematically in Fig. [2ja). 
This is essentially a spin echo experiment: it pulses are 
used to remove any net precession due to static magnetic 
fields, while rectifying the sinusoidal A-C phase accumu- 
lation. The disk is first set spinning. As the diamond 
passes point A, the NV center is optically pumped into 
the m = state |0) using the 532nm laser (Bloch sphere 
i of Fig[2](a)). Next, a ir/2 Rabi pulse creates a coherent 
superposition of |0) and |1) (ii). As the disk continues to 
rotate, a n Rabi pulse is applied every time the diamond 
passes points A or B, which inverts the phase between 
the |0) and |1) sublevels (iii,iv). After the diamond has 
completed n full rotations (v), a second 7r/2 pulse is ap- 
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Figure 2: (a) Pulse scheme (pulse times not to scale) and 
rotating-frame Bloch sphere state evolution for a single rota- 
tion of the disk, with n rad of accumulated A-C phase. Runs 
with multiple rotations would have further 7r pulses whenever 
the diamond passes points A and B. (b) Example fluorescence 
signal as a function of electric field strength, where the max- 
imum A-C phase is 10 rad. j/-axis units are arbitrary. Dots 
represent example experimental data points. 



plied, which converts the relative phase acquired during 
rotation into a population difference between the sub- 
levels (vi). A pulse from the 532nm laser then excites 
the NV center and a fluorescence measurement is taken. 
Note that we can ignore the motion of the diamond for 
the duration of each of the pulses (~0.5 pis or less): given 
reasonable figures for the disk's rotation frequency and 
radius (see below), the diamond is stationary to within 
a few tens of microns. The total phase accumulated $ is 
given by 



(12) 



An ensemble of measurements, made over time, will 
not be perfectly coherent. For NV centers, the dominant 
decoherence process is dephasing (divergence in the phase 
evolution of different members of the ensemble) rather 
than spin relaxation. Dephasing manifests itself as an 
exponential decay of the fluorescence signal, towards a 
value corresponding to equal populations of |0) and |1), 
as the evolution time of the experiment is increased. 

The fastest dephasing, due to inhomogeneous magnetic 
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fields, gives a so-called T 2 * dephasing time on the order of 
/xs [28|| . However, these magnetic fields fluctuate slowly 
compared with other time scales, which means that the 
sequence of ir pulses eliminates any net effect [l5| . The 
next fastest dephasing is the T2 homogeneous broadening 
due fast-fluctuating magnetic fields caused by the nuclear 
spin of 13 C atoms. T2 times of up to 1.8 ms have been 
recorded in recent experiments [171 ]. 

If the pulse sequence is conducted for £ ranging from 
zero to £0 while keeping the number of rotations per run 
n and the evolution time t r fixed, then the fluorescence 
count, plotted as a function of £ , will look like Fig. [2^b). 
To minimize error in measuring the maximum A-C phase 
from the fluorescence signal, the gradient of the curve 
should be a maximum at £q ■ This can be achieved by ap- 
propriately shifting the phase of the final n/2 microwave 
pulse. In the example shown, the pulse's phase lags by 
2.15 rad behind the earlier pulses. To measure a very 
small phase ($>ac "C 1 rad), the lag should be tt/2. 

Van Oort has shown that, in addition to the well- 
known excited state Stark shift [231 . the NV ground state 
experiences a linear Stark shift [301 ] . with Hamiltonian 

ffstark - -h£R 2E ( e ~ 34e |-l)(l| +e 3i9 |l)(-l|) , (13) 

where the constant R2E is approximately 20 HzV _1 cm 
and 9 is the angle between £ and one of the crystal's 
three planes of symmetry containing the N-V axis. The 
energy eigenstates then rotate at triple the frequency of 
the diamond's rotation. By the adiabatic theorem, pro- 
vided the diamond's rotation frequency is much smaller 
than the energy separation between the m — ±1 states, 
the Stark shift will simply shift the energy of these states 
by a constant amount, Y? R\ E £ 2 j (B z g[iB) ■ The 7r-pulses 
will therefore eliminate any net effect of the ground state 
Stark shift. 

We envisage realistic experimental parameters / = 4 
kHz, r = 1 cm, £ = 30 kV/mm 31], and T 2 = 1.8 ms. 
These parameters would give a total A-C phase of ~ 17 
rad, which is two orders of magnitude larger than previ- 
ous neutron and atom beam experiments 2, 11 1 . Exper- 
iments measuring the resistance of mesoscopic rings have 
demonstrated A-C-like phases of this order though 
in these cases the phase comes from precession about a 
continuously changing axis, as the electric field is perpen- 
dicular to the plane of motion of the magnetic dipole, and 
the Rashba spin-orbit coupling arises from a many-body 
interaction. 

Two intrinsic sources of error are shot noise, due to the 
finite fraction of photons detected, and spin projection 
noise, due to the probabilistic nature of quantum mea- 
surements. These both follow Poisson statistics, yielding 
a sensitivity r\ in the measurement of Qac, 



C (sa 0.05 for typical experiments) is a factor which takes 
into account collection efficiency and the nonzero fluores- 
cence of |1) [IH and T is the combined time of the indi- 
vidual runs. We have taken the duration t r of each run to 
be T2, to optimize the relative uncertainty A^ac/'&ac 
given the signal's exponential decay in t r jTi- 

A single NV center could measure an A-C phase to 
within a radian over 100 hours (77 w 700 rad Hz -1 / 2 ). 
Alternatively, a spatial ensemble of N individual centers 
would allow a high- precision measurement, improving the 
sensitivity by y/N. NV densities as high as 10 n mni" 3 
could be used without admitting spin-spin decoherence 
[IH which, given a crystal size of (1 mm) 3 , would give 
very high sensitivity: 77 w 2 mrad Hz -1 / 2 . 

We have shown that an A-C phase shift in diamond 
can be observed using current experimental techniques. 
The phase shift can be very large and can be measured to 
great accuracy. Moreover, as a solid state system whose 
path can be precisely controlled, the NV center can be 
used to probe deeper aspects of topological phases, such 
as the non-Abelian dynamics of the A-C phase in three 
dimensions. 
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